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Abstract 

This paper is devoted to study the energy content of colliding plane 
waves using approximate Noether symmetries. For this purpose, we 
use approximate Lie symmetry method of Lagrangian for differential 
equations. We formulate the first-order perturbed Lagrangian for col- 
liding plane electromagnetic and gravitational waves. It is shown that 
in both cases, there does not exist non-trivial first-order approximate 
symmetry generator. 
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1 Introduction 

For the description of a physical system, the conserved quantities like en- 
ergy, linear and angular momentum etc. are of great importance. Being a 
non-conserved quantity, energy and its localization is the most prominent 
and fundamental issue in general relativity. For gravitational waves (the 
ripples produced by accelerating mass in the fabric of the spacetime which 
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propagates through the space at the speed of light pQ), this problem is of 
particular interest. The gravitational waves, by definition, have zero stress- 
energy tensor. Thus the existence of these waves was questioned. However, 
the theory of general relativity indicates the existence of gravitational waves 
as solutions of the Einstein's field equations [2]. In fact this problem arises 
because energy is not well-defined in general relativity. The problem for 
plane gravitational waves was resolved by Ehlers and Kundt [3J, Pirani [I] 
and Weber and Wheeler [5] by considering a sphere of test particles in the 
path of the waves. They showed that these particles acquired a constant 
momentum from the waves. Qadir and Sharif |6j presented an operational 
procedure, embodying the same principle, to show that gravitational waves 
impart a momentum. 

To find a satisfactory solution of this challenging problem, enormous ef- 
forts have been made. The pseudo-tensors which are firstly introduced by 
Einstein provide a mean to determine global energy-momentum conservation. 
This attracted many others physicists to develop different prescriptions e.g. 
[Tj-JTS]. Being a combination of T\ and a pseudo-tensor t b a (as it explicates 
the energy and momentum density of the gravitational field), none of these 
non-tensorial complexes are proved to be unambiguous. 

Recently, an alternative approach, i.e., concept of approximate symmetry, 
is used to define energy of gravitational waves. There have been a number of 
definitions of approximate symmetry and different methods [HJ [T5] are avail- 
able to find these symmetries. Using this concept, a measure of the extent of 
the break-down of time-translational symmetry was presented [16]- [18]. This 
yields the so-called an almost symmetric space and almost Killing vector (the 
vector field corresponding to almost symmetric space) [19]. However, each 
of these attempts have their own drawbacks. 

Noether symmetries also known as symmetries of the Lagrangian, pro- 
vide a systematic treatment for the solution of DE arising in many practical 
problems [20]. Noether and Lie symmetries have many important applica- 
tions like the linearization of the non-linear equation, the reduction of the 
order of ordinary differential equations (ODEs) as well as the number of in- 
dependent variables of partial differential equations (PDEs) etc. The double 
reduction of DEs and the existence of conservation laws or first integrals 
(Noether invariants) are other interesting features of the Noether symme- 
tries [21]. The more symmetries an equation possesses, the more easier will 
be its integration. The field equations that provide a basis for general theory 
of relativity (GR) are highly non-linear coupled PDEs. Finding an exact 
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solution to these field equations is one of the major problems in GR. The 
Noether symmetries approach has been proved to be fruitful in this regard 
due to their wide range of applications in the field of DEs [221 123]. It is 
worth mentioning here that most of physically interesting solutions to the 
field equations, available in literature, possess some kind of symmetry [23]. 

The Killing vectors (KVs) form a subalgebra of Noether symmetries which 
further form a subalgebra of Lie point symmetries. By investigating Noether 
symmetries and the corresponding conservation laws for some static models, 
it was shown that Noether symmetries yield some non-trivial conservation 
laws that are different from KVs [22]. The perturbed Lagrangian yields the 
approximate symmetries and the approximate conserved quantities as con- 
stants of motion. By contracting the first order non-trivial approximate sym- 
metry with the momentum 4- vector, one can obtain energy non-conservation 
due to time variation [25]. Since Noether theorem leads to the fact that each 
continuous symmetry generator of the Lagrangian corresponds to a conserved 
quantity. Therefore, in order to define energy, it is very useful to utilize the 
concept of time symmetry (due to the fact that energy conservation comes 
from the time translational invariance). 

In going from Minkowski spacetime (flat) to non-flat spacetimes, some 
of the conservation laws are lost. Using the connection between spacetime 
isometries and its related differential equations (geodesies) [261 EZ] for differ- 
ent spacetimes, the lost conservation laws are recovered and the correspond- 
ing energy re-scaling factors are obtained [2S]-[2I]. We have used this pro- 
cedure to discuss energy in the stringy charged black hole solutions [32J and 
Bardeen model [33J. In a recent paper [25], energy contents of pp and cylin- 
drical gravitational waves have been discussed by using the same procedure. 
In this paper, we extend this procedure to colliding plane electromagnetic 
and gravitational waves to check the existence of conserved quantity. 

The plan of this paper is as follows. In the next section, we review the 
mathematical framework for exact and approximate symmetry methods for 
the Lagrangian of differential equations (DEs). Sections 3 and 4 contain ap- 
proximate Noether symmetries of colliding plane electromagnetic and gravi- 
tational waves respectively. The last section provides summary and results. 
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2 Mathematical Formulation 



Noether's theorem was proved in 1915 and published in 1918 by Emmy 
Noether. This theorem describes the connections between continuous sym- 
metries of a physical system and their corresponding conserved quantities 
[34] . By the analysis of various transformations, Noether's theorem provides 
amazing insights into any general theory. This theorem provides informa- 
tion regarding the conservation laws. If the results of a physical experiment 
are independent of any change in position (homogeneity in space) and time, 
then the formulated Lagrangian is symmetric under continuous spacetime 
translations and this theorem leads to laws of energy and linear momentum 
conservation [551 ES]- Similarly, if the physical experiment is independent of 
any change in the measurement of angle, then the physical system is rota- 
tionally symmetric and this theorem implies the conservation law of angular 
momentum. 

Noether's theorem has a wide range of applications in theoretical physics 
and calculus of variation. Dissipative system is one of the examples of the 
systems which cannot be modeled with Lagrangian. Therefore, Noether's 
theorem is inapplicable to these systems and consequently their continuous 
symmetries provide no conservation laws. The original form of this theorem 
is valid for those Lagrangian which contain only the first-order derivatives. 
In many practical problems, DEs arise with a small term called the pertur- 
bation parameter which indicates a small error or correction. The purpose 
for introducing this small term is to examine the given DE in some certain 
limit. The Lie point symmetries of such DEs (perturbed) are of great im- 
portance. The application of Noether symmetry analysis to Lagrangian of a 
system provides those symmetries which directly yield our desired conserved 
quantities. In that case, one needs to use only the first-order prolongation 
of the symmetry generator. The procedure for calculating the symmetries of 
the Lagrangian J37J EH] is given as follows. 

We take a vector field X as 



XW = X + ( V b s + r, b a x a - ^ s x b - C, a x a x b ) — , (a, b = 0, 1, 2, 3). (1) 




and its first-order prolongation is defined as 
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Suppose a second-order ODE, i.e., Euler-Lagrange equation is given by 

x a = g(s,x a ,x a ). (2) 

The vector field X is said to be Noether point symmetry of the Lagrangian 
L(s, x a , x a ) (corresponding to Eq.([2])), if there exists a function A(s, x a ) such 
that the following condition is satisfied 

X. [1] L+(D s OL = D S A. (3) 

Here D s is the total derivative operator given by 

and A is a gauge function. 

The significance of Noether symmetries is given by the following theorem. 

Theorem: Suppose L(s, x a , x a ) be the Lagrangian corresponding to second- 
order ODE given by Eq.(j2]) and X is its corresponding Noether point sym- 
metry. Then the first integral of motion associated with the Noether point 
symmetry X is defined (39] as 

I = tL + W-x°0—-A. (5) 

Now we define approximate Noether symmetries for DEs. The first-order 
perturbed Lagrangian corresponding to the first-order perturbed DE 

E = E + eEi + 0(e 2 ) 

is given by 

L(s, x a , x a } e) = L (s, x a , x a ) + eL^s, x a , x a ) + 0(e 2 ) (6) 

such that the functional j v Lds is invariant under the one-parameter group 
of transformations with approximate Lie symmetry generator 

X = X + eXt + 0{e 2 ) 

and 

A = A + eAi 
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is the first-order perturbed gauge function. Also, here 



X b = 6>— + ; a = 0,l,2,3, 6 = 0,1. 

Then the conditions for calculating exact and first-order symmetries are given 
as follows 

X^Lq + p.£o)A> = D S A , (7) 
XS 1] L + Xj, 1 ^! + (Ds&Lo + (A^o)^ = D S A U (8) 

where X and Xi are the exact and first-order parts, respectively, of the sym- 
metry generator X. The perturbed Lagrangian always admits a symmetry 
eXo called the trivial symmetry. If, for the symmetry generator X to exist 
with X 7^ and Xi ^ &X , where k is any arbitrary constant, then X is 
called the non-trivial symmetry. The corresponding exact and the first-order 
approximate parts of the first integral are given by 

BT 

h = CoL + (Vo-i a Co)g^-A , (9) 

BT BT 
h = ^L 1 +^ 1 L + (^-x%)-^ + ( V a l -x^ 1 )-^-A 1 . (10) 

The first integral of motion I = I + el\ is said to be unstable if I = 0; 
otherwise it is stable. 

For further details, one can see the literature 139 



2.1 Exact and Approximate Symmetries of the Collid- 
ing Plane Electromagnetic Waves 

Bell and Szekeres, in 1974, presented an exact solution of field equations 
due to the collision and the subsequent interaction of two electromagnetic 
waves (12] . They considered the case of two-step electromagnetic waves which 
collide on the flat background region. The spacetime is defined in terms of 
four regions [13]. Under the coordinate transformations u = v = 
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the spacetimes in these four regions can be written as 

Region(I) (t < z, t < -z) : 

dt 2 dz 2 
ds 2 = — dx 2 — dy 2 — , (11) 

Region(II) (t > z, t < -z) : 

ds 2 = ^-- cos 2 a ( t -^)(dx 2 + dy 2 )-^-, (12) 
Region(III) (t < z, t > -z) : 

ds 2 = ^ - cos 2 b^^dx 2 + dy 2 ) - (13) 

Region(IV) (t > z, t > —z) : 

dt 2 dz 2 
ds = — cos (At- Bz)dx 2 -cos 2 {Az- Bt)dy 2 , (14) 

where A = and B = The solution has a special property of being 
conformally flat, i.e., all components of the Weyl tensor vanish. 

In order to discuss Noether symmetries for colliding plane electromagnetic 
waves, we define Lagrangian for the spacetime given by Eq.(ll4p as follows 



L = ^ - cos 2 (At - Bz)x 2 - cos 2 (Az - Bt)y 2 - y. (15) 

To evaluate exact Noether symmetries of this Lagrangian, we substitute val- 
ues of the Lagrangian, the first-order prolongation of the symmetry genera- 
tor (given by Eq.flTJ) and the total derivative operator (given by Eq.(J3])) in 
Eq.([7|). Then we compare coefficients of the coordinate derivatives and their 
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products on both sides. In this way, we obtain the following set of DEs 



O = 0, e.<B = 0, i v = 0, C = 0, A iS = 0, (16) 

77° = A ttl -2cos 2 (At - Bz)r,] s = A tX , (17) 

-2 cos 2 (Az - Bt) v 2 s = A >y , -n* = A jZ , (18) 

77° - 2 cos 2 (At - Bz)rj] t = 0, (19) 

Tj° y - 2 cos 2 (^ - Bt)rj 2 = 0, 77° - ^ = 0, (20) 

cos 2 (At - Bz)^ + cos 2 (Az - Bt)r, 2 x = 0, (21) 

V % + 2 cos 2 (At - Bz)i)\ = 0, 277° - O = 0, (22) 

77J + 2 cos 2 (^ - Bt) V 2 z = 0, 2^ - C,s = 0, (23) 

2ri] x - £ a - 2(4 - B) tan(At - B;z)(?7 + 77 s ) = 0, (24) 

2 rfy -is~2(A- B) tan(Az - Bt) (77° + t/ 3 ) = 0. (25) 

This is a system of 19 DEs which we have to solve for the five unknowns 
£, 77 , 77 1 , r] 2 , 77 s using back and forth substitutions (where each of these is 
a function of five variables s, t, x, y, z). The solution (exact symmetry 
generators) to the above set of DEs can be written as 

S dt S dz ' 1 dx 1 2 dy ^ ^ 

It has been mentioned earlier that the electromagnetic waves collide 
with each other on the flat background region (i.e., Minkowski spacetime 
in Cartesian coordinates) and is represented by non-static spacetime (given 
by Eq. (jl4p ). Therefore we define it as a perturbation over some static space- 
time, i.e., spacetime representing the interaction region (region (IV)) as a 
perturbation over the flat Minkowski spacetime. The exact static spacetime 
is given by 

dt 2 dz 2 
ds 2 = ^L- dx 2 -dy 2 - -^- (27) 

and the corresponding Lagrangian can be written as 

L= t ~x 2 -y 2 -^. (28) 

We evaluate Noether symmetries of this Lagrangian by using Eq. (J7J) . Af- 
ter substituting all the corresponding values, we compare coefficients of the 
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coordinate derivatives and their products and obtain the following set of 
determining equations 



it = 0, i x = 0, i y = 0, £* = (), A >s = 0, (29) 

v° a = Au -2»?i = A»> -v 3 s = A^ (30) 

-2^ = ^, <,-2^ = 0, (31) 

7)1-2^ = 0, rj^- V 3 t = 0, r,l + V l = 0, (32) 

^. + 2^ = 0, 2^-^ = 0, (33) 

^ + 2^ = 0, 2^-^ = 0, (34) 

2^-^ = 0, 2^-^ = 0. (35) 

The simultaneous solution for this set of DEs is given by 

6 = -y- + c i s + c 2, (36) 

?7 = t( " C ° S ^~ °^ + sc 3 + 2xc 5 + 2c 7 y + c 12 z + c 13 , (37) 

i x(cos + Ci) 

77 = h sc 4 + te s + cay + ci 4 2 + ci 5 , (38) 

^ 2 = y{cos + ci) + s ^ + ^ + c ^ _ ^ + c ^ ^ 3Q ^ 

^ 3 = z(cqs + ci) _ gCg + fei2 _ 2xCu _ 2yci 6 + en, (40) 

X 2 Z 2 



X, y, *) = c (j - y - — - — ) +tc 3 - 2xc 4 - 2yc 6 + C 9 Z + C10, 

(41) 

where all c's are arbitrary constants of integration. Clearly this is a 17 di- 
mensional algebra in which many symmetries are present which provide no 
conservation laws. In this algebra, 10 symmetries correspond to the gen- 
erators which form the Poincare group and 7 provide the other significant 
generators which have been discussed in detail [30]. These symmetries can 
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be re-arranged to the following form 



st d sx d sy d sz d 
~ ~2di + Ydx + Y dy + ~2 ~b~z 

Yi = ^, Y 2 = 4, Y 3 = 24 + tA, 
s at at ax 

9 ^ T d __ „ 9 9 
Y 4 = 5 — , Y 5 = — , Y 6 = 2y— + t— , 
ay ox at ay 

d d d d 

ax ay az a^ 

d d d 

Yl ° = ^ + Yll = ^' 

9 9 a 
Y12 = 2x— , Y13 = — , 

ax az ax 

Y14 = z— 2y— , Y 15 = — , Y 16 = s— . 

ay oz ay ox 

In order to calculate conserved quantity by using the non-trivial first- 
order Noether symmetry (if it exists), we define the perturbed spacetime as 
follows 

dt 2 dz 2 
ds 2 = —- [1 + ecos 2 (A - Bz)]dx 2 - [1 + e cos 2 (Az - Bt)]dy 2 - — . (42) 

(Basically it means that we are taking the spacetime of region (IV) as a 
perturbation over the spacetime of region (I)). The corresponding first-order 
perturbed Lagrangian is given by 

i 2 z 2 
L = — - [1 + e cos 2 {At - Bz)]x 2 - [1 + ecos 2 (As - Bt)]y 2 - — . (43) 

The first-order approximate symmetries of this Lagrangian is found by using 
Eq.flH]). After substituting the Lagrangian and other corresponding values, 
we make comparison of the coefficients of coordinate derivatives and their 
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product terms. Consequently we obtain the following set of DEs 

O = o, O = o, iy = o, = 0, 4« = °> ( 44 ) 

7?° = A t , -2^ - 2cos 2 (At - Bz)[^°- + c 4 ] = A >x , (45) 

-2T? 2 - cos 2 (Az - Bt)[^ + ce] = A„, = A z , (46) 

?7° - 277J - 2c 5 cos 2 (At - Bz) = 0, (47) 

77° - 2?7 2 t - 2c 7 cos 2 (Az - St) = 0, 77° - r/J = 0, (48) 

^ + V 2 X + c 8 [cos 2 (At - Bz) - cos 2 (Az - Bt)} = 0, (49) 

r}% + 2^ + 2ci4 cos 2 (At - Bz) = 0, 2rf t - £ a = 0, (50) 

rjl + 2 V 2 Z + 2c 16 cos 2 (Az - Bt) = 0, 2 V 3 Z - £, s = 0, (51) 

277^ - i s - 2{A - B) cos(At - Bz) sm{At - &)[!^±fll 
^ z(c s + d) + ^ _ ^ + 2 ^ _ ^ + 2 ^ _ ^ 

+c i2 (t + z) + cis + en] + (c s + c x ) cos 2 (At - Bz) = 0, (52) 
2rf, l- O - 2(A - B) cos(Az - Bt) sm(Az - Bt) f( c ° s + c ^ 



H h s(c 3 - c 9 ) + 2x(c 5 - c u ) + 2y(c 7 - c i6 ) 

+ci 2 (t + z) + ci 3 + cu] - (c s + ci) cos 2 (Az - Bt) = = 0. (53) 

In this system of DEs, 14 out of 17 constants corresponding to exact sym- 
metry generators are present. We solve this system of DEs by back and 
forth substitutions and check whether there exist some non-trivial parts of 
the symmetry generators or not. Since all the constants corresponding to the 
exact symmetry generators given by Eqs. fl36l) - fj40l disappear, therefore the 
system of DEs given in Eqs. (|44p - (|53|) becomes homogeneous. Thus it yields 
symmetries of the static spacetime only. 

2.2 Exact and Approximate Symmetries of the Collid- 
ing Plane Gravitational Waves 

The spacetime representing the collision of two plane gravitational waves in 
Rosen form is given by 

ds 2 = 2e- M dudv - e- u {e v dx 2 + e~ v dy 2 ), 
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where U, V and M are functions of the null coordinates u and v. For the 
collision of such waves, one has the following division of the spacetime into 
four regions. The spacetimes in these four regions are given by [45J: 

Region(I) (u < 0, v < 0) : 

ds 2 = 2dudv - dx 2 - dy 2 , (54) 
Region(II) (u > 0, v < 0) : 

ds 2 = 2(1 ± u^-^dudv - (1 ± uy- a dx 2 - (1 ± u) 1+a dy 2 , (55) 
Region(III) (u < 0, u > 0) : 

ds 2 = 2(1 ± v) L ^ 21 dudv - (1 ± v) 1 ^^ 2 - (1 ± v) 1+a dy 2 , (56) 
Region(IV) (u > 0, u > 0) : 

ds 2 = 2(1 ± u ± v^-^dudv - (1 ± u ± v)^ a rfx 2 

-(1 ± M ±v) 1+a ^ 2 , (57) 

where a is any arbitrary parameter. If one chooses a = and the negative 
signs only, then this solution is quite similar to the solution obtained by 
Szekeres and Khan and Penrose [16]. Using the suitable coordinate transfor- 
mations u = v = ^ [44J and taking only the negative signs, we have 
the spacetime in region (IV) as follows 

ds 2 = 2(1 - t)^ - ^) - (1 - t)(dx 2 + dy 2 ). (58) 
To evaluate exact symmetries of the Lagrangian for this spacetime 

L = 2(1 — t)^Cj ~ j) ~ (1 - t)(x 2 + y 2 ) (59) 

we use Eq.flTJ) and substitute the values of first-order prolongation operator, 
the total derivative operator and the above defined Lagrangian in it. Then 
we compare the coefficients of the derivatives and their products, we obtain 
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the following set of DEs 

O = o, O = o, e,» = 0, c, = o, As = o, (60) 

77° = At, -2(1 - *hi = A tX , (61) 

-2(i-of. = A». -^7f== = A- (62) 

77° -2(l-t) 3 / 2 77j = 0, (63) 

^0 _ 2(1 _ t)S/2^J = o, 7/1-^ = 0, (64) 

^ + ^ = 0, + 2(1 - = 0, (65) 

2(l-*)[27/J-£,,]+T7° = 0, v 3 y + 2(l-tf 2 V 2 z = 0, (66) 

2(1 - *)(2»fi - e,.) - »7° = 0, (67) 

(l-t)[2 V l x -U-V° = 0, (l-t)[2 V l-U-V° = 0. (68) 

The simultaneous solution of this system of DEs turns out to be 

v ° = 0, v 1 = b 6 y + b 7 , rj 2 = -b 6 x + b 8 , rf = b b , £ 2 = 6 2 , A = 64 ,(69) 

where all 6's are constants, which show that the gauge function is constant. 
These can be re-arranged to the following form 

d d d d d 

A = b , £2 = c , Yi = — , Y 2 = y- x— , Y 3 = — , Y 4 = — . 

az ox ay ox ay 



In order to discuss the corresponding conserved quantity, we take Eq. (158]) 
as a perturbation over background static spacetime in a similar pattern as we 
have done for colliding electromagnetic waves. The background static space- 
time is the Minkowski spacetime in Cartesian coordinates given by Eq.( l271) 
and its Noether symmetries are given by Eqs. fl36|) - fj40|) . T evaluate the ap- 
proximate Noether symmetries, the first-order perturbed spacetime for the 
collision of gravitational waves is given by 

ds 2 = [1 + e(l - t)^){^ - ^) - [1 + 6(1 - t))(dx 2 + dy 2 ). 
The corresponding Lagrangian is formulated as 

L = [1 + 6(1 - f)^](£ - y) - [1 + e(l - t)](x* + y*). (70) 
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Using this Lagrangian and the corresponding values in Eq.flS]), we obtain the 
following system of DEs 

It = 0, £,* = 0, i y = 0, e i2 = 0, A iS = 0, (71) 

V° s + (f + e>)-±= = A* (72) 

-2^-2(l-t)(^ + c 4 ) = A x , (73) 

-2^-2(l-t)(^ + c 6 ) = A„ (74) 

-^--^=(^-c 9 ) = A tZ , (75) 

K - H] + ^f= t c ^ 1 ~ ^ ~ *) 3/2 ) = °' ^ 
K - H\ + jr^pi^ 1 - ^ - 1 ) 3/2 ) = °. ( 77 ) 

V° z -V 3 t = 0,vl + V 2 x = 0, (78) 
fe + 2^]--|^=(l-(l-t) 3 / 2 ) = 0, (79) 

iH - U + 2(1 _ 1 t) 3/2 ^( C g + C + SC 3 + 2XC 5 + 2c ty 

+c 12 z + c 13 ] = 0, (80) 
^ + 2^-^(l-(l-t) 3/2 ) = 0, (81) 

(2< - y + 2 (i-t)3/2 ^ c ° g + Ci ) + SC3 + 2xC5 + 2c7?/ 

+ci 2 2; + ci 3 ] = 0, (82) 

[2^ - is] - (^(cqs + ci) + sc 3 + 2xc 5 + 2c 7 ?/ 

+c 12 ^ + c 13 ) = 0, (83) 

[2^ - i s ] - (-(c s + ci) + sc 3 + 2xc 5 + 2c 7 y 

+c 12 z + ci 3 ) = 0. (84) 

In this system of 19 DEs, 12 out of 17 constants, corresponding to exact 
symmetry generators defined in Eqs. (!36|) - (|4*0|) . appear. When we solve this 
system of DEs simultaneously then in the final solution all these constants 
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disappear. Therefore it makes the system of equations homogeneous and no 
non-trivial symmetry generator exists. 

3 Summary and Discussion 

In this paper, we have discussed energy contents for the colliding plane waves. 
For this purpose, we have applied the slightly broken or approximate Lie sym- 
metry methods for Lagrangian to colliding plane electromagnetic and gravita- 
tional waves. Since the contraction of non-trivial approximate Noether sym- 
metry (corresponding to time) with the momentum 4-vector can be used to 
define energy imparted to test particles [22] ■ Therefore, we have defined the 
first-order perturbed Lagrangian for both these solutions and check whether 
there exists some non-trivial symmetry. 

We have obtained a system of 19 DEs both for the colliding plane elec- 
tromagnetic as well as colliding plane gravitational waves. For the colliding 
plane electromagnetic waves, 14 out of 17 constants whereas for colliding 
plane gravitational waves, 12 out of 17 constants corresponding to exact sym- 
metry generators defined in the Eqs. fl36|) - fj4*0|) appear. When we solve these 
systems of DEs simultaneously, all these constants disappear in both cases. 
Thus, the resulting systems of DEs becomes homogeneous and are the same 
as that for the exact case (i.e., Minkowski spacetime in Cartesian coordi- 
nates). Consequently only the previous exact symmetry generators of static 
spacetime which correspond to 17 dimensional Lie algebra for Minkowski 
sapcetime in Cartesian coordinates are obtained. These symmetry gener- 
ators provide the conservation laws for energy, linear momentum, angular 
and spin angular momentum and hence there does not exist any non-trivial 
symmetry. 

Thus we conclude that for these first-order perturbed Lagrangian for col- 
liding plane electromagnetic and gravitational waves, there does not exist 
the conserved quantity. It is quite similar to the case of pp and cylindrical 
gravitational waves [25] where no non-trivial symmetry generator exists and 
one cannot determine the conserved quantity. It is mentioned here that one 
cannot apply these methods to the case of electromagnetic waves as there 
is no self-interaction. However, it may give significant result about the self- 
damping or enhancement of the gravitational waves if one applies this method 
to the system of geodesic equations for these spacetimes. 
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